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Abstract 

We have developed a compact Compton gamma-ray imager with a large held of view and a low channel-count that 
is capable of quickly localizing gamma-ray sources in the few hundred keV - several MeV range. The two detector 
planes (scatter and absorber) employ bars of Nal(Tl) read out by photomultiplier tubes (PMTs) located at each end. 
The long-range imaging performance has been tested from 392 keV to 1274 keV. An angular resolution measure of 
2.72° ± 0.06° and an efficiency of (1.79 ± 0.04) x ICE 3 at 662 keV is obtained. A 137 Cs (662 keV) source equivalent 
to a 10 mCi source 40 m away can be located in 60 seconds with an uncertainty of about a degree. No significant 
degradation in imaging performance is observed for source angles up to 40° off axis. 

Keywords: Compton Imaging, Compton Telescope, Nal, Scintillator, Security, Radiation Detection, Pulse-Height 
Sharing 


Acknowledgements 

This is NRCan/ESS Contribution 20140132. 

1. Introduction 

A portable Compton gamma-imaging device is a helpful tool in the prevention or clean-up of a radiological or 
nuclear incident. First responders require portable radiation detectors operating in the 0-3 MeV range that can be 
deployed in the held. To meet their needs, we have developed a Compton imager that is low-cost, ruggedize-able, and 
adaptable to suit a number of different surveying platforms (truck, helicopter, etc.). The detector is highly sensitive, 
(capable of localizing a 10 mCi 137 Cs source at 40 m distance with degree-level precision in 60 s), making it suitable 
for real-life situations involving radiation sources that can be weak, distant or shielded. 

Some highly-compact Compton-imaging designs based on scintillator flTj and semiconductor technologies (CZT |2] 
|3])) have been developed for commercial use. In the future, development in crystal growth technology may make these 
viable options for affordable larger-scale detectors. However, at present, these small ultra-portable detectors have lim¬ 
ited sensitivity. In parallel, extremely large highly-sensitive imagers have been developed BE 012). However, 
these do not meet the requirements of agencies with limited funds, nor are they suitable for mounting onto multiple 
platforms. Our detector aims to operate in the middle range, by being both sensitive and portable (can be loaded into 
and out of a vehicle), while remaining affordable. A more expensive design based on scintillator with SiPMs has been 
developed J8). The design presented in this paper employs long bars of scintillator (Nal(Tl)) with photomultiplier 
tube (PMT) readout, a combination well-proven for application in robust low-cost detectors operating in the field. 


*now at the Geological Survey of Canada. 

Email address: audrey. macleod@nrcan-mean. gc. ca (A.M.L. MacLeod) 


Preprint submitted to Preprint submitted to Nuclear Instruments and Methods in Physics Research A 


June 18, 2015 






Figure 1: Compton imager schematic depicting the energy and position of the scattered electron, E\, jci, recorded in the scatter plane, and the 
energy and position of the scattered gamma ray, £2 and * 2 , recorded in the absorber plane. The Compton cone, of opening angle 6c, is calculated 
using the energy deposits and localizes the possible source locations. 


2. Compton Imaging 

In a Compton gamma imager, an incoming gamma ray of energy E y Compton interacts in a scatter detector, 
depositing energy E\. The outgoing scattered gamma ray deposits its energy £2 in an absorber detector. Neglecting 
Doppler broadening, the scattering angle between the initial and final-state gamma rays, 9 C , can be determined from 
the two energy deposits, via 

, 1 1 

cos e c = 1 + m 0 c~(— - —), (1) 

where E y - E\ + £2 and moc 2 is the electron rest energy. As illustrated in Figure [I] the position of the source lies 
somewhere on a cone of opening angle 8c, where the cone axis is on the line between the two energy deposits. The 
location of the source can be reconstructed from the intersection of several cones. 

3. Bar Design 

We have designed, constructed and tested a Compton imager made of long bars of Nal(Tl) read out by PMTs fixed 
to the bar ends. Pulse-height sharing between the PMT signals is used to determine the position of the interaction. 
The use of long bars of Nal(Tl) offers an inexpensive design solution, as Nal(Tl) is relatively low cost, and the bar 
format requires few readout channels. 

Long scintillation bars have previously been employed in Compton imagers where the energy deposits are large [0, 
however to employ this technique at lower energies, an optimization of the attenuation length of the bars is needed. 
Studies have been done to look at the effect of tuning the attenuation length of scintillator in order to optimize its 
position reconstruction mmmm, which holds well at high incident gamma-ray energies where one is not photo¬ 
statistics limited. However, operating at low energies (less than 150 keV in the scatter layer for incoming 662 keV 
gamma rays), requires tuning the bar attenuation length to optimize the overall performance of the imager. To this end, 
a highly detailed Monte-Carlo simulation using optical transport was developed to determine the optimal attenuation 
length of the scintillator bars. A prescription was then developed to manufacture these bars for the detector. 

4. The Detector 

Previous work has shown that for the scatter plane, the optimal thickness for one Compton scatter of 662 keV 
gamma rays in Nal(Tl) is ~ 20 mm 03]. With this thickness, approximately 20% of 662 keV gamma rays Compton 
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Figure 2: Prototype consisting of ten scatter modules (front) and seven absorber modules (back) with a scatter-absorber spacing of 425 mm. 



Figure 3: A single absorber detector module. The scintillation light from a bar of encapsulated Nal(Tl) is read out by PMTs placed at the ends of 
the bar. The bar and PMTs are enclosed in a thin aluminum shell and packed with foam for protection. 


scatter once and exit. For the absorber plane, a thickness of 40 mm is chosen, a compromise among increasing the 
absorption fraction, maintaining comparable lateral and longitudinal position resolution, and matching the bar cross 
section with commercially available PMTs. 

A photograph of the full-scale prototype is shown in Figure [2] The modular scatter and absorber layers are 
425 mm apart (center to center) and mounted on a frame constructed of 1/2” PVC tubing. The scatter plane con¬ 
sists of ten 16 x 16 x 200 mm 3 and 20 x 20 x 200 mm 3 bars of Nal(Tl) and the absorber plane consists of seven 
40 x 40 x 400 mm 3 bars of Nal(Tl). 

An overview of the scintillator dimensions and the PMT sizes, QEs and model numbers of all the components are 
provided in Table [T] In the scatterer, we used both Super bialkali (SBA) PMTs, with QE ~ 35%, and Ultra bialkali 
(UBA) PMTs, with QE ~ 43%. In the absorber, we used standard PMTs. The scintillator was purchased from Scionix, 
Saint Gobain and Hilger Crystal^] while the PMTs were from Hamamatsif] 

Each Nal(Tl) bar is encapsulated in aluminum with a quartz window located at each end. The bars are read out by 


2 Scionix: www.scionix.nl; Saint Gobain: www.saint-gobain.com; Hilger Crystals: www.hilger-crystals.co.uk. 
3 www.hamamatsu.com 


Plane 

Scintillator 

PMT 

number 

size (mm 3 ) 

model 

size (mm) 

QE 

Scatter 

5 

20 x 20 x 200 

R8900U-100 

23.5 

SBA (35%) 

4 

16 x 16x200 

R8900U-100 

23.5 

SBA (35%) 

1 

16 x 16x200 

R7600U-200 

18.0 

UBA (43%) 

Absorber 

7 

40 x 40 x 400 

R6236 

54.0 

23% 


Table 1: Summary of the detector components making up the scatter and absorber plane. Only four of the seven absorber bars were included in the 
detector for assessing the performance for reasons discussed in Section[6] 
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PMTs glued to the windows. The entire setup is positioned inside a 1/16” thickness aluminum container (see Figure|3]i 
and packed with foam to protect the assemblies from mechanical shock. The PMT high voltage (HV) is supplied by 
a CAEhfjSY 2527 Universal Multichannel Power Supply System. 

The data acquisition uses a combination of VME, CAMAC and NIM standards. The PMT pulses are amplified 
by a factor of ten (Phillips ScientificJ^] 776) and split into two sets. The first set goes into a discriminator (Phillips 
Scientific 705) and participates in a hardware coincidence effectively requiring one or more end-to-end coincidences 
in the scatter layer together with the same pattern in the absorber layer (Phillips Scientific 754 Logic Units). The 
second set of amplified PMT signals is fed directly into a 48-channel 10-bit VME VF48 digitizer M and integrated 
on-board when a coincidence is present. The resulting charge values are then transferred to a computer via a USB- 
VME bridge (CAEN V1718). Immediately, a software cut is made to include only events with exactly one scatter bar 
hit and one absorber hit, where a hit denotes that both the left and right PMTs of a given bar are above threshold. The 
energy deposit is reconstructed in a given bar by summing the signals of the left PMT (SI) and right PMT (S2). The 
position is reconstructed using the ratio of these two signals, S1/S2. 

5. Development and Performance of Scintillator Bars 

5.1. Motivation 

The photon propagation in a scintillator is determined by its physical and optical properties as well as its surface 
finish, coating and wrapping. The key to our design is tuning the attenuation length of the scintillator bars in order to 
improve the position resolution while maintaining as good an energy resolution as possible. The attenuation can be 
tuned by modifying the surface finish, e.g. polishing or degrading the smoothness mmm. This section provides a 
detailed explanation of how the energy and position reconstruction depend on the attenuation. The optical absorption 
coefficient, a, multiplied by the length of the bar, L, is used here as a measure of the attenuation. 


Position Resolution 

"x . 
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Figure 4: Schematic of the position reconstruction using a bar of length L indicating S i and S 2 , the signals from the left and right PMTs, and x, 
the displacement of the interaction from the bar center. 


5.2. Position Resolution vs. Attenuation Length 

For an interaction at the longitudinal position x in a bar of length L, the signal S 1 from the left PMT can be written 
as 


and the signal from the right PMT as 


Si = 


N 0 


Pe 


-a(i+x) 


( 2 ) 


S 2 


— Pe~ a ^' 
2 


■x) 


(3) 


where No is the number of photons generated and P is the probability of detecting each photon. The probability, P. 
incorporates all sources of photon loss that are independent of the position of interaction, such as the probability of 


1 www.caen.it 

4 http://www.phillipsscientific.com 
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generating a photon, transmission at the scintillator ends, and PMT QE. L and x are the scintillator length and the 
displacement of the interaction from the bar center, respectively (see Figure[4|. 

Propagating the photostatistical uncertainties of S \ and S 2 , yields an expression for the position resolution fTTI : 

cr Xsm = — 7 ==e^ V cos Hax). (4) 

aVMrP 

Equation [4] highlights what is known to be true empirically. The position resolution depends on the initial number of 
photons produced, No, the probability of detecting each photon, P, and the absorption coefficient, a. The minimum 
of the position resolution along a bar of length L occurs when aL = 2.9 CD. aL for a polished bar of Nal is much 
less (aL ~ 0.3 for a 20 x 20 x 200 mm 3 bar). Thus, increasing a should improve the position resolution. However, 
increasing a also has the undesirable effect of increasing the position dependence of the position resolution. 

5.3. Energy Resolution vs. Attenuation Length 

The expression \jS \ x .S' 2 is proportional to No and independent of x and can thus be used to reconstruct the energy 
without knowing the position where deposition occurred. A derivation similar to the position resolution calculation can 
be used to obtain an expression for the relative uncertainty of the energy due to statistical fluctuations by propagating 
the statistical uncertainties of S \ and S 2 - 


°~ £s ““ = e 4 yfcosh(ax). (5) 

E VW 5 

Equation [5] shows that increasing the number of detected photons, either by increasing No or P, improves the 
energy resolution whereas increasing a worsens the energy resolution. Increasing a also has the undesirable effect of 
increasing the position dependence of the energy resolution. The situation is therefore complicated. Decreasing the 
attenuation length (increasing aL ) both improves the position resolution and worsens the energy resolution. Thus, to 
determine the ideal attenuation length, an optimization is undertaken with simulations. 

5.3.1. Relating the Energy Resolution to the Position Resolution 

The total energy resolution, cr E , is equal to a combination of cr £stal , the energy resolution due to statistical fluctua¬ 
tions, and (Tem, the intrinsic energy resolution, 


2 _ 2 2 
O' E — ^"Estat + O' Elnt • 


(6) 

By combining Equations [4] and [5j one obtains an expression relating the statistical component of the position 
resolution to the statistical component of the energy resolution, 


_ 1 tr £ sta 

0"Xs M [ --— 

a E 


(7) 


Equation[7]is then combined with Equation[6]to obtain an expression for the position resolution, <r x , as a function 
of the measured total energy resolution, (^f), 


* crxs, “ “ l V(?)M^r)- (8) 

By assuming a fixed value of 6% for the intrinsic energy uncertainty, we use Equation[8]to obtain an estimate of the 
position resolution indirectly from the energy resolution measurement at energies for which the position resolution 
is not directly measured. Note that the intrinsic energy resolution depends on the quality of the crystal and is not 
measured for our scintillator. Thus, the assumption of a fixed 6% intrinsic energy uncertainty contributes to a moderate 
error in the estimation of the position uncertainty at higher energies (~ 5% error at 662 keV), with the error decreasing 
at low energies where the uncertainties from photostatistics dominate. 
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5.4. Optimal aLfrom Simulations 

GEANT4 1731 simulations were used to determine the optimal attenuation length of the bars. 

First, a detailed optical-photon transport was carried out to determine the energy and position resolution for 
Nal(Tl) bars as a function of attenuation length. Different attenuation lengths were modeled by making part of the 
bar surface absorbent. The surface treatment is referred to as polished front-painted in GEANT4 and is specularly 
reflective with 96% reflectivity. The bars simulated were of sizes 20 x 20 x 200 mm 3 and 40 x 40 x 400 mm 3 and 
had 0%, 25%, 50%, 75% and 100% of one side surface-treated, exploring aL values in the range of ~ 0.3 to 2.5. A 
monoenergetic 662 keV gamma-ray line source was directed towards the bar at multiple positions along its length, and 
the reconstructed energy and position resolution values were used to extract resolution functions as input for further 
modeling. 

Simulations were then carried out for a complete detector made up of ten 20 mm-width scatter bars and seven 
40 mm-width absorber bars, with a vertical spacing between the bars of 65 mm for the scatterer and 90 mm for 
the absorber. The separation between the scatter and absorber plane was fixed at 425 mm. The performance of 25 
detectors corresponding to the different combinations (5 x 5) of scatterer and absorber were compared. For each 
detector combination, a simulation equivalent to a three-hour exposure to a 10 mCi 137 Cs source located on the 
detector symmetry axis 40 m away was carried out. The energy and position of the interactions recorded in GEANT4 
were smeared using the energy and position resolution functions derived from the detailed optical simulation. Only 
events with full gamma-ray energy deposit shared between exactly one scatter bar and exactly one absorber bar were 
selected. Furthermore, events exhibiting the energy signature of a back-scatter event (incoming gamma scattering ~ 
180° backwards from absorber to scatterer) were rejected. 

The width of the angular resolution measure (ARM) for the selected events was used as a figure of merit for 
detector performance. ARM = 0 C - 0 geom , where 0 (: is the Compton angle calculated from the energy deposits 
and 0 g eom is the observed scattering angle based on the known source position and the reconstructed hit positions. 
Figure [5] (left) shows the standard deviation, cr ARM , as a function of scatter bar aL, with the absorber aL fixed at 1.01 
(25% of one side surface-treated). A significant improvement (~ 26%) in cr arm is observed, decreasing from 3.1° 
to 2.8° with increasing surface treatment. Figure [5] (right) shows ct A rm as a function of absorber bar aL, with the 
scatter aL fixed at 1.14 (50% of one side surface-treated). Once again, <x ARM decreases with increasing aL, leveling 
off at aL ~ 1. The ARM for the remaining possible combinations of scatter and absorber aL demonstrated similar 
behaviour. Since decreasing the attenuation length (or increasing a) also increases the non-linearity of the energy and 
position resolution along the bar, an aL of 1, where the improvement in <x ARM levels off, was deemed optimal. 
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Figure 5: Simulation. Left: as a function of scatter bar aL, for absorber aL fixed at 1.01 (25% of one side surface-treated). Right: o"arm as 

a function of absorber bar aL, for scatter aL fixed at 1.14 (50% of one side surface-treated). Note that the y axis is zero-suppressed and the error 
bars are smaller than the symbols in both the left and right figures. 
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5.5. Determination of the Surface Treatment 



Figure 6: Plot of the ratio, (S\) / (S 2 ) for a 20 X 20 X 200 mm 3 bar at 662 keV, as a function of position for a polished bar (blue stars); four sides 
degraded with 220 grit sandpaper (green triangles); two sides degraded with 220 grit and the remaining two sides with 120 grit (red squares); and 
four sides degraded with 120 grit (black circles). The data sets can be described by exponential functions. Note that the data have been scaled such 
that the middle point is located at (0,1). 


Surface Treatment 

a L 

A E 

E 

(% FWHM) 

cr x 

(mm) 

Polished 

0.23 

7.3 

18 

220 

0.90 

8.4 

7 

220/120 

0.92 

9.8 

8 

120 

1.14 

9.5 

6 


Table 2: Scionix data for aL measured for a 20 mm bar with successive levels of sanding and the energy resolution measured at the bar center at 
662 keV. The position resolution values presented here are estimated from the energy resolution and a using the technique described in Section|5.3.1| 


A technique to modify the bar surface was developed with the manufacturer Scionix to achieve scintillator bars 
with attenuation lengths of aL ~ 1. Due to the hygroscopic nature of Nal(Tl), any alteration of the scintillator surface 
must take place in a low-humidity environment. Thus, the procedure to modify the bar surface was undertaken in a 
dry room at the manufacturer’s site and carried out by the manufacturer. 

To tune the attenuation length, the bar surface was degraded by sanding with various grit sizes in the longitudinal 
direction. Four different surfaces were compared: polished; four sides degraded with 220 grit sandpaper; two sides 
degraded with 220 grit and the remaining two sides with 120 grit; and four sides degraded with 120 grit. Figure [6] 
shows the ratio of the mean values of the signals at the two ends, (S 1 ) / (S 2 ), versus the position of a collimated 
662 keV energy source along the bar, for these four surface treatments. As the sandpaper grit number decreases 
(becomes more abrasive), the attenuation of the bar increases. The attenuation length was obtained by applying an 
exponential fit to the data. As shown in Table [2] aL increases from 0.23 for an untreated bar to 1.14 for 120 grit 
applied to all four sides. 

Energy resolutions measured with the 662 keV source at the center of the treated bars are also presented in Table[2j 
along with position resolutions calculated using the technique described in Section 5.3.1 (Note that the ratio values 
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used to determine aL and the energy resolutions were provided by Scionix without uncertainties.) 

We find, as expected, a dramatic improvement in position resolution as aL increases to close to 1 with only 
minor improvements with further degradation of the bar surface. The accompanying worsening of energy resolution 
is comparatively minor. The second most attenuating surface treatment (220 grit to two sides and 120 grit to the two 
other sides) was chosen for the 20 mm bars, yielding aL = 0.92, which is close to the targeted value of 1. 

The same procedure was then repeated for the 16 mm and 40 mm bars to obtain aL values ~ 1. For the 16 mm 
bar, the same combination of 220 grit and 120 grit was chosen, while for the 40 mm bar, a surface treatment of 120 
grit was chosen. To test the reproducibility of the technique, bars were then ordered from two other manufacturers, 
Hilger and Saint-Gobain, who employed the same prescription for the surface treatment. 


5.6. Performance of the Bars 

In total, 18 bars were tested from three manufacturers, Scionix, Hilger, and Saint-Gobain. The bars came in three 
sizes: 16 mm width (five bars), 20 mm width (six bars), and 40 mm width (seven bars). Each set had the same surface 
treatment applied to it. For in-house testing, the bars were paired with PMTs having areas that matched their cross 
section. The 16 mm and 20 mm bars were paired with R8900U-100 PMTs (QE = 35%) while the 40 mm bars were 
paired with R6236 PMTs (QE = 23%). The bars were tested in their aluminium containers described in Section[4] 

We made a detailed measurement of the performance of the bars at a single energy, 662 keV. We parametrized the 
position as a function of the left and right PMT signals to reconstruct the position of the interactions. A 10 /rCi 137 Cs source 
was placed behind a cylindrical collimator having a thickness of 70 mm and an aperture of 5.5 mm. The collimator 
was attached to a motorized platform (Arrick Robotic^Model XY-30). Four scans of each bar were made, two in each 
direction, with 5 mm steps. The left and right PMT outputs (S \ and S f) were obtained by digitizing the PMT signals 
with a 1 GS/s, 8-bit digitizer (Acqiri^]DC270) and computing the area of each pulse. The digitizer was triggered by 
coincidence of the left and right PMT signals. At each step, 10 000 coincidence events were recorded. 

To obtain the position response, Gaussians were fit to the photopeaks in S i and S 2 and the photopeak pulse-height 
means, (S 1 ) and ( S 2 ), were measured. Figure [7] shows the signal ratio, (S 1 ) / (S 2 ), versus distance along the bar, for 
an individual 20 mm bar from each of the three manufacturers. 

To determine the position resolution, only events within 2.5 standard deviations of the mean of the photopeaks 
in S 1 , S 2 , and their sum, were considered. For these events, the ratio, (Si)/ (S 2 ), was calculated and converted to 
a position using the position versus signal characterization. The width of the resulting distribution of positions was 
adopted as the position resolution at 662 keV. 

Calibration of the bars at other energies was performed differently. Sources were placed only on the bar center. A 
Gaussian was fit to the summed signal to obtain a parametrization of the energy resolution and energy versus signal. 

To obtain a parametrization of the position resolutions for these energies, the approach described in Section 5.3.1 was 
used. 

The average aL, and the energy and position resolutions obtained at 122 keV and 662 keV for the 16 mm and 
20 mm scatter bars are shown in Tables [3] and [4] respectively. The uncertainties on the attenuation measurements 
(when N, the number of bars measured, was greater than one) were estimated from the spread between the different 
bars and their relatively small values (~ 5%) indicate that the application of the surface treatments is repeatable by 
a given manufacturer. The bars from the different manufacturers have similar attenuation. However, the measured 
energy and position resolutions differ significantly from one manufacturer to another, with the Hilger bars exhibiting 
the worst performance. 

The performance of the bars was also compared with the performance of a simulated bar of equivalent aL. These 
results are shown in Table[3]for the 16 mm bar, and similar results were observed for the other bar sizes. The measured 
energy and position resolution values at 662 keV were found to be in good agreement with the simulations. However, 
those for 122 keV were found to be slightly worse than predicted (<20%), with the Hilger bars deviating the most 
from the simulations. The Hilger bars exhibited significantly worse energy and position resolution values. For these 
bars, we observed far fewer photoelectrons, which could help explain why both the energy and position resolution 
values are worse. This demonstrates the limitations in the simulations to model the optical photon transport of the bars 


5 www.arrickrobotics.com 

6 www.agilent.com 
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of Nal(Tl) accurately, where the properties of the scintillator, scintillator surface, and PMT photocathode are difficult 
to model and also vary from one manufacturer to another. 

The average aL, and energy and position resolutions obtained at 662 keV for the 40 mm absorber bars are shown 
in Table [5] The Scionix bars are significantly more attenuating (aL = 1.10) in comparison to those from the other two 
manufacturers {aL ~ 0.8). 



Position (mm) Position (mm) Position (mm) 

Figure 7: Plot of the ratio, (Si)/ (S 2 ), for an individual 20 mm bar from three manufacturers (Saint-Gobain (left), and Hilger (middle) and Scionix 
(right)). To obtain aL, the central region of each set of points is fit to an exponential function extending from - 60 mm to + 60 mm. 


Manufacturer 

-Vbar 

aL 

^ 122 keV 

^ 662 keV 

cr x 122 keV 

cr x 662 keV 




% (FWHM) 

% (FWHM) 

(mm) 

(mm) 

Hilger 

2 

1.02 + 0.06 

17.5 ± 0.2 

8.1 ±0.3 

13.1 ±0.9* 

5.6 ± 0.4 




(13.3 + 0.3) 

(8.4 ± 0.6) 

(10.8 + 0.1) 

(5.2 ± 0.3) 

Scionix 

3 

1.18+0.01 

15.4 + 0.4 

8.2 ±0.1 

11.3 + 0.2* 

5.1 ±0.1 




(13.7 ± 0.4) 

(8.5 ± 0.7) 

(9.6 ±0.1) 

(4.9 ± 0.3) 


Table 3: aL, and the energy resolution and position resolution measured at both 662 keV and 122 keV for the 16 mm bars provided by Hilger and 
Scionix. A^g^Rrepresents the number of bars tested. The position resolution values followed by * are estimated from the energy resolution using 
the technique described in Section [5.3.1| The performance of a simulated bar of equivalent aL is presented in parentheses. 


Manufacturer 

Near 

aL 

122 keV 

% (FWHM) 

M 662 keV 

% (FWHM) 

cr x 122 keV 

(mm) 

cr x 662 keV 

(mm) 

Saint-Gobain 

1 

0.863 ± 0.001 

13.8 + 0.1 

8.0 ± 0.2 

12.1 ±0.1* 

5.7 ±0.1 

Hilger 

2 

0.999 ± 0.003 

18.4 + 0.6 

8.7 ± 0.2 

14.7 + 0.5* 

6.2 ±0.3 

Scionix 

3 

1.03 + 0.02 

15.6 + 0.5 

7.7 ± 0.3 

11.9 + 0.6* 

6.0 ±0.1 


Table 4: aL, and the energy resolution and position resolution measured at 662 keV and 122 keV for the 20 mm bars. The position resolution 
values followed by * are estimated from the energy resolution using the technique described in Section |5.3.1| 
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Manufacturer 

Near 

aL 

662 keV 

% (FWHM) 

cr v 662 keV 

(mm) 

Hilger 

2 

0.76 ± 0.03 

12.3+0.9 

19.6 + 2.7 

Saint-Gobain 

3 

0.78 ± 0.02 

8.2 ± 0.2 

10.9 ± 0.9 

Scionix 

2 

1.07 + 0.08 

7.9 ±0.1 

8.5 ±0.1 


Table 5: aL, the energy resolution and position resolution measured at 662 keV for the 40 mm bars. 


6. Bar Failures 

Three of the seven absorber bars experienced a sudden change in performance. The attenuation length decreased, 
accompanied by an increase in the number of photoelectrons detected at one end of the bar and a decrease at the 
other end. Both the Scionix and Saint Gobain bars had this issue, but not the Hilger bars. This is currently under 
investigation and there are several possible hypotheses, such as contamination of the Nal crystal, problems with the 
glue on the glass-scintillator interface, failure of the plastic seal that surrounds the quartz windows at the bar ends, and 
cracking due to the thermal expansion of the crystal, glass and aluminum. Integrating the PMT into the encapsulation 
is a possible development that would improve ruggedness and may mitigate this problem. 

7. Detector Performance 


Energy 

(keV) 

Source 

Source Strength 
s“' (xlO 6 ) 

392 

113 Sn 

3.2 ±0.1 

662 

137 Cs 

22.5 ± 0.5 

1274 

22 Na 

22.3 ± 0.5 


Table 6: List of the radioactive sources imaged along with the gamma-ray energy and source strength. While 22 Na provides two peaks, 511 keV 
and 1274 keV, only the imaging performance of the 1274 keV peak was assessed. 


Due to the failure of three absorber bars, only four absorber bars, located centrally in the absorber plane, were 
included in the detector performance studies. 

Long-range imaging tests were performed in a large room (16mx 17 m) with a false floor at the National Research 
Council (NRC) in Ottawa, Canada. A list of the sources used is provided in Table[6] The gamma-ray emitting sources 
were placed 9 m away from the detector, at different angles [{(f> = 0°, 9 = 0, 10, 20, 30°), (<p = 180°, 9 = 10, 20°) 
and (<p = 90°, 9 = 10, 20°)]. The origin of the coordinate system used lies at the center of the scatter plane. The 
x-axis lies along the length of the bars, the y-axis points upwards, and the z-axis points out of the front-face of the 
detector away from the absorber. The azimuthal angle, <p, lies in the x-y plane and 9 is the polar angle with respect to 
the z-axis. 

7.1. Event Selection 

A description of the trigger selection is provided in Section[4] Several criteria were used to select events off-line. 
First, events with exactly one scatter bar and exactly one absorber bar above threshold were selected. Second, a ~ 2-cr 
energy cut around the photopeak in the total energy distribution was implemented. Third, to reject back-scatter events, 
where a gamma ray is scattered in the absorber first, and then absorbed in the scatterer second, a ~ 2-cr exclusion cut 
was applied to the scatter-energy back-scatter peak. Fourth, a minimum threshold of 30 ke V was required in the scatter 
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Source 

Energy 

(keV) 

Scatterer 

(keV) 

Sum 

(keV) 

113 Sn 

392 

30 - 100 

340 - 440 

137 Cs 

662 

30-150 

600 - 720 

22 Na 

1272 

30- 180 or 280-650 

1150- 1350 


Table 7: The event selection used includes a back-scatter rejection cut and a photopeak selection cut. 



Figure 8: Scatter energy plotted against summed energy for 137 Cs (662 keV). The vertical lines indicate the cuts performed on the summed energy. 
The horizontal lines indicate the selection performed on the scatter energy. The events that fall into the rectangular region formed by the overlap 
of the four lines were selected. Note that events from the back-scatter peak, which is visible at a total energy of 0.66 MeV and a scatter energy of 
0.19 MeV, were thereby excluded. 


detector. A summary of the energy selection used for the three sources is provided in Table [7] A visual representation 
of the cuts used for a 137 Cs source is shown in Figure]^ 

Figure |9] shows the summed energy distributions (scatter + absorber) for a I37 Cs source before and after the 
selection criteria are applied, overlaid with the estimated contribution from naturally occurring radioactive material 
(NORM). The NORM contribution was determined by taking a run of equal duration with no source present. Energy 
distributions for the 22 Na and 113 Sn sources are shown in Figure[l0] The varying fractional contribution from NORM 
to the total energy distribution is due to the different source strengths. There also exists a significant number of 
accidental coincidence events with NORM. 

We define the detector efficiency as N/Nj 0Vd \, where N is the number of events that satisfy the event selection 
criteria and Njotai is the total number of gamma rays passing through the active area of the front face of the detector 
(active area = 40 000 mm 2 ). To determine the efficiency, we relaxed the photopeak cut and in its place used a fit 
to isolate events in the photopeak. All of the other selection criteria were applied. To account for the number of 
background events, a Gaussian and a third order polynomial is fit to the summed energy distribution, see Figure [TT| 
The correctly reconstructed forward-going scatters fall in the Gaussian part of the distribution. The efficiency of the 
detector for a 662 keV source placed off-axis (6 = 10°) was measured to be (1.79 ± 0.04) x 10 3 . The efficiencies 
at 392 keV and 1274 keV are provided in Table[8] 
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Figure 9: 1,7 Cs (662 keV) energy distribution. White - all data. Grey Shaded - estimated contribution from NORM. Red Lined - selected events. 




Figure 10: Energy distributions for —Na (left) and ll3 Sn (right) from 2 hours of exposure. White - all data. Shaded Grey - estimated contribution 
from NORM. 
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Figure 11: The summed energy distribution for 137 Cs (662 keV) fit to a Gaussian and a third order polynomial. 


Energy 

Efficiency 

(keV) 

(xl(T 3 ) 

392 

2.00 ± 0.09 

662 

1.79 + 0.04 

1274 

1.10 + 0.03 


Table 8: Efficiency for three sources positioned at 6 = 10°, (f) = 0°. The uncertainties reflect both the statistical fluctuations and the uncertainty in 
the source emission rate. 
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7.1.1. Angular Resolution Measure 




Figure 12: ARM distribution for a 137 Cs source located on-axis at 6 = 0° for selected events comparing measurement (left) and Monte Carlo (right). 
The distribution is fit to the sum of a Gaussian and a third-order polynomial. The widths of the Gaussian fits are shown. 

Figure [l2] (left) shows the ARM distribution for the selected events of a 137 Cs source located on-axis. The ARM 
distribution has been fit to the sum of a Gaussian distribution and a third-order polynomial. The width of the Gaussian 
fit yields cr arm of 2.72° ± 0.06°. 

The experimental ARM was compared with GEANT4 simulation. The detector configuration matched the imager 
that was tested, with ten scatter bars and four absorber bars. The simulation made use of realistic energy and position 
resolutions based on experimental data (see Section |5.6| l. The recorded energies and longitudinal coordinates (x) 
of the interactions were smeared by Gaussian probability distributions with widths determined by the energy and 
position resolutions measured for each bar. Both the z and y-coordinates of the interactions were assumed to be at 
the bar centers. The simulated ARM distribution is shown in Figure 12 (right), and the width of its Gaussian fit is 
cr arm = 2.66° ± 0.03°. The good agreement with experimental data confirms our understanding of the detector 
performance. 


7.2. Back-Projection Image 

A back-projection image is a 2-D distribution formed by back-projecting Compton cones into angle space. Fig¬ 
ure 13 shows a back-projection image from ~ 20 000 Compton cones (two hours of data) collected with the 137 Cs source 


located 9 m away and 20° off axis. The back-projection image correctly indicates the position of the source. 
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X (degrees) 


Figure 13: A back-projection image of selected events for a l37 Cs source. The x and y-axis are presented in vertical and horizontal degrees. The 
true source position is indicated by cross-hairs. The broad excess at the origin is caused by NORM. 


7.3. x 2 Minimization Algorithm 

An iterative ^-minimization procedure using the MINUIT package m was employed to determine the vector 
s(0,<b) which best represents the direction from the origin to the source. For N events, a x 1 function is constructed 
with two fit parameters, 6 and <p , representing the polar and azimuthal angles of the source direction vector which are 
to be determined: 


£[ACA(W,f 

jri °ACA/ 

where ACA (6, cj)) is the angle of closest approach of the direction vector s to cone i, and <x ACA , is the uncertainty of 
this angle. The expression for ACA is given by 

ACA = arccos(s(0, (/>) ■ r,) - df , (10) 

where f, is the unit-vector axis for cone i, and 6 C is the Compton scattering angle of the ;th cone. The uncertainty 
of ACA is calculated on an event-by-event basis from the cone-axis uncertainty and the Compton opening angle 
uncertainty (using Equation 1), taking into account the uncertainties on the energy deposits and their positions. 

To extract the source direction from a given sample of N events, three iterations are performed, with the direction 
determined in each iteration passed on as a starting seed direction for the subsequent iteration. The highest bin of the 
coarse-binned back-projection image is used as a seed in the first iteration. For a detailed explanation of the algorithm 
see ref jTTI . 
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7.4. Imaging Performance at 662 keV 

To characterize the imaging performance of the detector, the 17 Cs source rates at a 9 m distance were converted 
into an equivalent 10 mCi source at 40 m with 100% branching ratic]^] To quantify the ability to reconstruct 6 and cf>, 
the events were divided into groups or trials of 60 s of data, and the ^-minimization algorithm was applied to each 
group. To determine the accuracy, the means of the resulting distributions of fitted 8 and <p were used as a measure 
of the average reconstructed source location. To determine the angular precision, the 2-D standard deviation of the 
reconstructed direction was used. Table[9]summarizes (0). ((/>) and the angular precision obtained from sixty 60-second 
trials. The sources are all reconstructed to within about a degree of the known source position. The small discrepancy 
observed in the reconstructed angles versus the known source position can be explained by the fact that the quoted 
uncertainties on the reconstructed angles are statistical in nature and do not include a number of factors such as the 
effect of NORM, trigger threshold variations, and the depth of interaction in the bars. The angular precision obtained 
is on the order of a degree. 
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Figure 14: Angular precision for an equivalent 10 mCi 662 keV source located 40 m away as a function of integration time for five source locations 
(0=0°; 6 = 0°, 10°, 20°, 30°, and 40°). At all angles, the angular precision reaches one degree in about a minute. Note that the data points are 
correlated and each point includes the entire data set. 


8 The relative source strengths are taken into account, as is the ^ dependence and the attenuation by the air. 
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Run # 

True Source Position 

Reconstructed Source Position 


e 0 

<t>o 

(0) 

(<P) 

Angular Precision 


(°) 

(°) 

(°) 

(°) 

n 

1 

0.1 

95.2 

0.3 ± 0.2 

- 

1.1 ±0.1 

2 

10.0 

179.1 

10.2 + 0.1 

180.0 ± 0.4 

1.0 ±0.1 

3 

20.0 

179.6 

21.4 + 0.1 

180.3 ± 0.2 

1.1 ±0.1 

4 

9.9 

0.7 

9.6 ± 0.2 

-0.9 ± 0.5 

1.2 ±0.1 

5 

20.0 

0.3 

19.7 + 0.1 

-0.9 ± 0.2 

1.0 ±0.1 

6 

30.0 

0.2 

30.5 + 0.1 

-0.8 ± 0.2 

1.1 ±0.1 

7 

40.0 

0.2 

39.6 + 0.1 

-0.7 ± 0.2 

1.1 ±0.1 

8 

9.9 

90.7 

9.3 ± 0.2 

91.1+0.9 

1.3 ±0.1 

9 

20.1 

89.3 

20.6 ± 0.2 

89.2 ± 0.4 

1.4 ±0.1 


Table 9: Summary of the fit values obtained at 662 keV: (6), (0) and the angular precision for 60-second trials. Note that the source position for 
Run #1 is very close to the z-axis, where 0 is undefined, thus the 0 value obtained for this run is not meaningful. 


Figure [14] shows the angular precision at 662 keV as a function of time for several source positions. The detector 
demonstrates a wide field-of-view with good performance for sources located up to 40 ° off-axis. 

7.5. Performance at Different Energies 

The 1274 keV and 392 keV peaks of a 22 Na and 113 Sn source located at 6 - 10°, <!> - 0° were also successfully 
imaged. The results were again scaled such that the rates of gamma rays incident on the detector are equivalent to that 
from a 10 mCi source at 40 m with 100% branching ratio, allowing for direct comparison with the 662 keV data. Due 
to the significant contribution from NORM to the 113 Sn data, it was necessary to set the starting seed position used by 
the imaging algorithm to the known location of the source. For 60-second runs, we obtained an angular precision of 
1.0° + 0.1° and 1.6° ± 0.2° for 22 Na and 113 Sn, respectively. 

8. Conclusion 

We have developed a Compton-scatter gamma-ray imager composed of bars of Nal(Tl). The key to the low-cost 
design is the use of pulse-height sharing to reconstruct events, made possible by the successful development of a 
surface treatment that was applied to the bars. The technique was shown to be robust and the bars were incorporated 
into a fully-functional imager. The detector demonstrates good imaging performance over a wide range of energies 
(392 keV, 662 keV and 1274 keV) and over a wide field of view. An angular precision of about one degree has been 
achieved in one minute for a 10 mCi 662 keV gamma-ray source with 100% branching ratio, located 40 m away. 

We are currently investigating the use of different scintillator materials (Csl) and more compact readout devices, 
leveraging the experience our group has in constructing pixellated Compton imagers based on CsI(Tl) read out with 
silicon photomultipliers mm. 
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